Broken axisymmetry phase of a spin-1 ferromagnetic Bose-Einstein condensate 
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A spin-1 ferromagnetic Bose-Einstein condensate subject to a certain magnetic field exhibits a 
broken-axisymmetry phase in which the magnetization tilts against the applied magnetic field due to 
the competition between ferromagnetism and linear and quadratic Zeeman effects. The Bogoliubov 
analysis shows that in this phase two Goldstone modes associated with U(l) and SO(2) symmetry 
breakings exist, in which phonons and magnons are coupled to restore the two broken symmetries. 
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I. INTRODUCTION 

Bose-Einstein condensates (BECs) with spin degrees 
of freedom have attracted growing attention since the 
first observation of the spin-1 23 Na BEC by the MIT 
group [1, 2]. In contrast to a magnetic trap, in which 
hyperfine-spin degrees of freedom are frozen, an opti- 
cal trap can confine atoms in all magnetic sublevels of 
spin, allowing study of the magnetic properties of BECs. 
A variety of experiments have been performed to date, 
on areas such as spin domains [3], interdomain tunnel- 
ing [4], and realization of a spin-2 23 Na BEC [5]. The 
spin-exchange dynamics of 87 Rb BECs have been inves- 
tigated experimentally by Schmaljohann et al. [6], Chang 
et al. [7], and Kuwamoto et al. [8]. 

Theoretical investigations of the spinor BEC have also 
been carried out extensively. Mean field theory (MFT) 
for a spin-1 BEC was formulated by Ho [9] and Ohmi and 
Machida [10]. The MFT of a spin-2 BEC was developed 
by Ciobanu et al. [11] and Ueda and Koashi [12]. Law 
et al. [13] developed a many-body theory of spin-1 an- 
tiferromagnetic BEC. Koashi and Ueda [14] and Ho and 
Yip [15] extended it to including the linear Zeeman effect 
and found that an antiferromagnetic BEC realize a frag- 
mented BEC for a weak magnetic field. The Bogoliubov 
analysis was carried out by Huang and Gou [16] and by 
Ueda [17] in the presence of the linear Zeeman effect. 
Their results agree with those obtained using a diagram- 
matic method by Szepfalusy and Szirmai [18]. In these 
studies, the Zeeman effects are restricted to those up to 
the linear order in the magnetic field. A unique feature 
of trapped atomic systems is that linear and quadratic 
Zeeman effects can be manipulated independently due to 
spin conservation. If we take the quadratic Zeeman term 
into account, the ground-state phase diagram becomes 
much richer as shown in Ref. [2]. In particular, under 
a certain range of linear and quadratic Zeeman effects, 
there is a special phase in which the magnetization tilts 
against the applied magnetic field. The investigation of 
some of the unique features of this phase is the primary 
purpose of our study. 

When a weak magnetic field is applied along the quan- 



tization axis, the m = 1 or — 1 state is favorable for 
a spin-1 87 Rb atom due to the linear Zeeman effect and 
the ferromagnetic interaction, where m refers to the mag- 
netic sublevel. On the other hand, the quadratic Zeeman 
effect raises the energy of the m = ±1 states relative 
to that of the m — state. As a consequence, if the 
quadratic Zeeman effect is sufficiently large, the spin vec- 
tor of the ferromagnetic ground state not only shrinks 
but also tilts against the direction of the magnetic field. 
Therefore, even if the Hamiltonian is axisymmetric with 
respect to the direction of the magnetic field, the ground 
state spontaneously breaks the axisymmetry. This phase, 
which we shall refer to as a broken-axisymmetry phase, 
was predicted in Ref. [2], but little attention has been 
paid to it from the viewpoint of axisymmetry breaking. 



In the present study, we investigate the Goldstone 
modes of this phase by studying its excitation spectrum. 
The BEC with ferromagnetic interactions has three 
phases: ferromagnetic, polar, and broken-axisymmetry 
phases. In the ferromagnetic and polar phases, only the 
U(l) (global phase) symmetry is broken, and the Gold- 
stone mode corresponds to a phonon. In the broken- 
axisymmetry phase, the SO (2) symmetry (axisymmetry) 
of the spin vector is broken in addition to the U(l) sym- 
metry. Because of the simultaneous breaking of these two 
continuous symmetries, the associated Goldstone modes 
are expected to involve both phonons and magnons. 



This paper is organized as follows. Section II reviews 
the mean- field ground state of a spin-1 BEC to make 
this paper self-contained. Section III uses the Bogoli- 
ubov theory to derive one gapful mode and two gap- 
less Goldstone modes. Section IV explores the implica- 
tions of the present study for other related studies, and 
Sec. V concludes this paper. Appendix A derives ana- 
lytic expressions for the broken-axisymmetry phase. Ap- 
pendix B discusses excitations in the ferromagnetic and 
polar phases for comparison with those in the broken- 
axisymmetry phase. 



2 



II. GROUND STATE WITH BROKEN 
AXISYMMETRY 

A. Formulation of the problem 

We consider a uniform system of N identical bosons 
with hyperfine spin 1 in which an external magnetic field 
is applied in the z direction. The Hamiltonian of the sys- 
tem is written as the sum of one-body part Hi and two- 
body interaction part Hu. The one-body part is given 

by 

Hi= I dr¥ m f-^V 2 -pm + qmAi> m , (1) 

m—— 1 ^ ' 

where subscripts m = +1,0,-1 denote the magnetic 
quantum numbers along the z axis, M is the mass of 
the atom, and p and q are the linear and quadratic Zee- 
man coefficients, respectively. In the case of spin-1 23 Na 
and 87 Rb atoms, q is positive. The two-body part, which 
is described by a contact-type s-wave interaction at ul- 
tralow temperature, takes the form 

Ku = \ 9f I ta*l,¥ m ,(m'-,ri\V F \m-,n)* m * n , 

F=0,2 J 

(2) 

where gF — 4-7rft 2 a_F /M with ao and a 2 being the s-wave 
scattering lengths in the singlet and quintuplet channels, 
respectively, and Vf projects a two-body state into that 
with total spin F. The absence of the projection onto 
the F = 1 channel is due to the Bose statistics. 

Because the system is uniform, it is convenient to 
expand the field operators in terms of plane waves as 
^ m = Q -1 / 2 J2q e tq ' r fflq,m) where f2 is the volume of the 
system and a qjm represents the annihilation operator of 
a boson with wavenumber q and magnetic sublevel m. 
Equations (1) and (2) are then rewritten as 

Hi =^2(ew- pm + qm 2 ) al m a Km , (3) 

k k 



4 P/|Ci|n 
l+1> 



q/|Ci|n 



|-1> 



FIG. 1: Ground-state phase diagram for a spin-1 ferromag- 
netic BEC. The dashed curves indicate second-order phase 
boundaries. In the figure, j + 1) and | — 1} represent the fer- 
romagnetic phase, and |0) represents the polar phase. The 
shaded region is the broken-axisymmetry phase, in which the 
magnetization tilts against the z axis. 



state is given by 

Wbec = ^ : (a ,m"o,m) : 
l 

+ (-pm + qm 2 ) 4, m a ,m - T^s, 

m— — 1 

(5) 

where s = (a 2 , — 2ao.iao.-i) /V3 is an annihilation op- 
erator for a singlet pair. In the MFT, we replace the 
operator do, m with a c- number CmVNo- Here, TVo is 
the number of condensed atoms and the order param- 
eters C m 's are complex variational parameters that are 
determined so as to minimize the energy functional un- 
der the constraint of normalization ^2 m \( m \ 2 = 1- For 
this purpose, we introduce a Lagrange multiplier \x and 
minimize (TL) — M-^o X) m ICm| 2 with respect to ( m . In 
the following, we denote the set of the order parameters 
as C = T (Ci, Co, C— l) ? where the superscript T stands for 
transpose. 



B. Ground states 



where e k = h 2 k 2 /{2M), c = (.go + 2.g 2 )/3, ci = 

(.92 - 5o)/3, p k = Eq,m a q,m«q+k,m, and f k = 

E q , m ,„ f m,™a q , m aq+k, n with f = (f x , f y , f z ) being the 
spin-1 matrices in vector notation. The symbol :: denotes 
the normal ordering of operators. The spin-spin interac- 
tion is ferromagnetic if ci < and antiferromagnetic if 
Ci > 0. It is known that the interaction between spin-1 
23 Na atoms is antiferromagnetic and that the interaction 
between spin-1 87 Rb atoms is ferromagnetic [19, 20]. 

Assuming that a macroscopic number of atoms occupy 
the k = state, we replace the relevant operators with 
c-numbers. The Hamiltonian for the BEC in the k = 



The ground-state phase diagram for a spin-1 ferromag- 
netic BEC is shown in Fig. 1 [2] . The phases are classified 
as follows: 

1. Ferromagnetic phase ( | + 1) and | — 1) in Fig. 1). 
The order parameter is given for p > by Cf = 
T (e ixi ,0,0) and for p < by Cf = T (0, 0^ e l *-i) 7 
where \m denotes an arbitrary phase of ( m , i.e., 

Cm \Cm |c ^ m . 

2. Polar phase (|0) in Fig. 1). The order parameter is 
given by Cp = T (0,e ix °,0). 
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FIG. 2: (Color online) (a) Transverse magnetization, i.e., 
magnetization perpendicular to the direction of the applied 
magnetic field {fx) = {Fx)/No as a function of linear 
and quadratic Zeeman coefficients. The transverse magne- 
tization is nonzero only in the broken-axisymmetry phase, 
(b) Schematic illustration of the spin vector in the broken- 
axisymmetry phase. 



We note that this ground state breaks the axisymme- 
try around the z axis despite the fact that the Hamilto- 
nian including the external magnetic field is axisymmet- 
ric. Thus, the ground state features spontaneous break- 
ing of axisymmetry or spontaneous breaking of the SO (2) 
symmetry. Such an axisymmetry breaking is due to the 
competition between the linear and quadratic Zeeman 
effects and the ferromagnetic interaction. The quadratic 
Zeeman effect decreases the z component of the spin vec- 
tor. However, a decrease in the length of the spin vector 
costs the ferromagnetic interaction energy. To reconcile 
the quadratic Zeeman effect with the ferromagnetic in- 
teraction, the spin vector tilts against the z axis. In fact, 
$ in Eq. (10) is a monotonically decreasing function of 
p and a monotonically increasing function of q, and the 
length of the spin vector (9) attains the maximum value 
of Nq for \c\\n — > oo. 



3. Broken-axisymmetry phase (shaded region in 
Fig. 1). The order parameter is given by (see Ap- 
pendix A for derivation) 



C±i = {q±p) 



IP 2 + 2 \ci\ nq - q 2 
8 |ci | nq 3 



Co 



V-P 2 ) (p 2 +2\c 1 \nq + g 2 
4 I ci I nq 3 



o*(X!+X-i)/2 



(6) 



In the broken-axisymmetry phase, the transverse magne- 
tization, which is perpendicular to the external magnetic 
field, 



{Fx) = ^(F x ) 2 + (F y y 



^q 2 ^^(p 2 + 2\c 1 \nq) 2 -q± 
- N r-j ; 2 ' ( 7 > 

2 |ci| nq 1 

is nonzero as shown in Fig. 2. (If we choose Co to be 
real and positive, we have (F x ) = (F±) coscfi and (F y ) = 
{Fx) sm<p, where (j) = Xi = — X-i-) The longitudinal 
magnetization, which is parallel to the external magnetic 
field, is given by 



(F z ) = N 



p(p 2 + 2\a\nq- q 2 



2\ Cl \nq 2 

The total magnetization is therefore given by 

\J Ac 2 n 2 q 2 — (p 2 — q 2 ) 2 



(8) 



(F)\ = \l (Fx) 2 + (F z ) 2 = N - 



2 \ci\ nq 



(9) 

The magnetization thus tilts against the applied mag- 
netic field with the polar angle 



-°- „f^^vV + 2Nn g + g n (iQ) 

V PVP 2 + 2\ci\nq- q 2 J 



III. BOGOLIUBOV EXCITATIONS AND 
GOLDSTONE MODES 



According to the Goldstone theorem [21], there exists 
a gapless excitation mode when a continuous symme- 
try is spontaneously broken. In the broken-axisymmetry 
phase, we have shown that the relevant continuous sym- 
metry is the SO (2) axisymmetry. Since we have treated 
the system in the MFT in which the global phase of the 
wave function is assumed to be arbitrarily chosen, the 
U(l) symmetry is also broken. Thus, the two continu- 
ous symmetries are simultaneously broken in this phase. 
In this section, we examine the corresponding Goldstone 
modes using the Bogoliubov theory. 



A. Basic theory 



We first formulate a number-conserving Bogoliubov 
theory [22] for a BEC with spin degrees of freedom [17]. 
The advantage of this number-conserving formulation is 
that we do not need to introduce the chemical potential 
as a Lagrange multiplier to adjust the particle number to 
a prescribed value. In this formulation, we replace arj, m 



With Cm (N - El, 



At 



1/2 



in Eqs. (3) and (4) 



and keep terms up to those of second order in dk^u,m an d 
&w 0to - We then obtain an effective Bogoliubov Hamil- 
tonian as 
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W e ff H ( £ k-pm + gm 2 -q(fz) -Cin + Ci7l|Co ~ 2&C-1 H «k,m«k,r, 

k^O m=-l 

+ cin</) • J] E /™,"4, m «k,„ + ^ Yl ( 2 ^k + ©k©-k + vivi^ 

k^O m,n k/O 

+ ^E(2i r I-> k + ^k-i r -k + ^I-i rt -k) +£o, 



k^O 



r 



(ii) 



where X> k = E ro C s k,m, ^k = E m , n /m,»C«k,«. and 
i?o represents a constant term. 

In general, for spin-/ atoms, we can express quasipar- 
ticle operators &k,<r's as linear combinations of the anni- 
hilation and creation operators of the original particles: 



B k = U(fc)A k + V(k)A*_ k . 



(12) 



Here XJ(k) and V(fc) are (2/ + 1) x (2/ + 1) real matrices 
and the bold letters represent sets of operators 

Ak = T («kj , ak,/-i , • • • , Ok,-/) , 

-^k = T (^k,/ ' ^k,/-l ) " " " ) ^k,-/) ' 

where crj is the label for each Bogoliubov mode. The 
quasiparticle operators (12) should satisfy the Bose com- 
mutation relations, 



0. 



<5k.k' 8 a ,a> 



(13) 



&k,<r, ^k',(T 

which lead to 

Y,[ U -M T U l ^{k)-V atl {k) T V^{k)\ = 5^, 



(14) 

J2[ U ^( k ) T V*^'(k)-V^(k) T U l , a ,(k)] = 0. (15) 

i 

We can rewrite Eqs. (14) and (15) in a matrix form, 

T (U + V)(U-V) =1. (16) 

Thus, U and V are not independent of each other. For 
later convenience, we write Bk as 



Bi 



1 



(u + v) (A k + A*_ k ) + (u - v) (A k - A*_ k 

(17) 



We seek the excitation spectrum E a and operators 6 kjCT 
such that the quasiparticles behave independently, i.e., 

K c ti = J2 E„blJ Ka + E vac , (18) 

k#0 <7={cr 1 ,<7 2 ,--- ,<T 2 / + i} 



where TL c s is given in Eq. (11), and E vac is the energy of 
the vacuum state for the quasiparticles. From Eq. (11), 
the Heisenberg equation of motion takes the form 

i/jAA k = M(fc)A k + N(fc)Al k , (19) 

where M(k) and N(fc) are real and symmetric (2/+ 1) x 
(2/ + 1) matrices. Using the quasiparticle Hamiltonian 
(18) and the commutation relations (13), on the other 
hand, we obtain 



ih— B k = E(fc)B k , 



(20) 



where E(fc) is the diagonal (2/ + 1) x (2/ + 1) matrix, 
whose diagonal elements correspond to the energies of the 
elementary excitations -Ek.o-j • Then substituting Eq. (19) 
into Eq. (20) and using Eq. (16), we obtain 

(M + N) (M-N) T (U + V) = T (U + V)E 2 . (21) 

Since E 2 is also a diagonal matrix, the Bogoliubov exci- 
tation spectrum can be found as the eigenvalues of the 
matrix 



G = (M + N) (M - N) 



(22) 



We note that G, which is the product of two Hcrmitian 
matrices, is not, in general, Hermitian. The present ap- 
proach has the advantage that we can reduce the dimen- 
sion of the eigenvalue equation from 2(2/+ 1) to (2/ + 1) 
and therefore the diagonalization is simplified. That is, 
instead of the diagonalization of the 2(2/+ 1) x 2(2/+ 1) 
matrix as 



M N 

-N -M 



E 

-E 



(23) 



a (2/ + 1) x (2/ + 1) matrix G is to be diagonalized. 



B. Low-lying modes in the broken-axisymmetry 
phase for k — > 

Without loss of generality, we may assume Cm to be 
real and positive. The excitation spectra in the ferro- 
magnetic and polar phases can be derived analytically as 



■5 



shown in Appendix B. The analytic solutions can also 
be obtained for the broken-axisymmetry phase. How- 
ever, since they are very complicated, we here derive the 
excitation spectrum for small k. 

The effective Hamiltonian (11) gives the coefficient ma- 
trices M and N of the Heisenberg equation of motion (19). 
Using the explicit form of ( m in Eq. (6), the matrix G 
can be written in the form, 

G = Go + 2(g 2 nG 1 - anG'Jek + Je£, (24) 

where / is the unit matrix and 

/ OiC-iCo 9iCiC-i ©iCiCo \ 

Go = OoC-iCo e CiC-i ©oCiCo , (25) 

\e_xC-iCo e_iCiC-i ©-iCiCo/ 
/ Ci 2 CiCo CiC-A 

Gi - CiCo Co 2 C-iCo , (26) 

VCiC-i C-iCo C 2 i / 

/C 2 C-i/Ci -2C-1C0 -2C1C-A 

Gi = -2C_iCo Co 2 + 2CiC-i -2C-iCo , (27) 

\-2CiC-i -2C-1C0 QCi/C-iJ 

with 

9 m = (3m 2 -2)\ Cl \n[p 2 + 2\ Cl \nq+(-l) m q 2 

-2mpq\( a /{q( m ). (28) 

I 



We first consider the limit of ek — > 0. The eigenvalues 
for G = Go can be obtained easily: one is 

E 2 gap = (3p 2 - 2 Cl nq - q 2 )(p 2 - 2c x nq + q 2 )/q 2 (29) 

and the other two are zero. Thus the system has two 
gapless excitation modes, which, as we will show later, 
arise from the U(l) and SO(2) symmetry breakings. We 
label the gapful mode as a, and the other two gapless 
modes as [3 and 7. 

The eigenvectors of Go are given by each row of the 
following matrix: 

/ A6i A9 A9_i \ 
U + V = (/ip + Mm)Ci MpCo (Mp - Mm)C-i , (30) 

\(>p + I'm)Cl ^pCo (^p-^m)C-l/ 

where A, /x p , /z m , v p , and v m are arbitrary parameters. 
Note that the second and the third rows are the lin- 
ear combinations of two basis vectors (Ci>Co)C-i) an d 
(Ci,0, — C-i)i both of which are eigenvectors with zero 
eigenvalue. It follows from Eq. (16) that the matrix U— V 
is given as the transposed inverse matrix of Eq. (30), 



C-1C0 

tj _ y _ 2_ I -C-iQo^p-^mHCoQ-i^p 



C-iei(^ P -^m)-Cie-iK+^m) 
j 

-C-lQl(Mp-Mm) + ClB-l(Alp+A'm) 
J 



C1C0 

A 

CiQo(' / p+' / m)-CoQi' / p 
J 

-Cieo(M P +Mm)+CoQiA'p 
J 



(31) 



where A = 2CiC_i0o 



CoC-i0i - CoCi©-i and J 



Low-lying modes in the broken-axisymmetry 
phase for small k 



In the limit of small ek, the five parameters A, /x p , 
v p , and v m can be determined by substitution of Eqs. (30) 
and (31) into the definitions of the Bogoliubov operators 

(17) and by comparison of the quasi-particlc Hamiltonian 

(18) with the effective Hamiltonian (11). However, we 
cannot perform this procedure using only the expression 
for ek = because the two Goldstone modes diverge in 
the limit of £k 0. Thus, it is necessary to find the ek 
dependence of the eigenenergies to find the properties of 
the low-lying excitations. From Eqs. (25)-(27), we obtain 



the eigenenergies up to the order of ek as 

'£ 2 =£ g 2 a P + 4^^)e k + 0(e2), 

< E 2 = A+e k + O (e k ) , 
E 2 =A_e k + 0(e 2 ), 



(32) 



where 



A± = g 2 n 




[2g 2 n - r)) 2 



8g 2 n(q - r]) v 2 
c\n(2>rj — 2q + 2c\n) 



with r] = (q 2 — p 2 )/q. In Fig. 3, we compare the ek de- 
pendences of the approximate eigenenergies (32) (dashed 
curves) with those of the numerically obtained exact en- 
ergies (solid curves). 

It is important to note that the two gapless excitations 
and E 1 in Eq. (32) share the same leading-order term 

1/2 

e k ' . Since the effective Hamiltonian H e s in Eq. (11) 

contains only the terms that arc proportional to ek, this 
1/2 

e k ' dependence must be canceled by the operators 6k, <r 
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0.2 0.4 0.6 0.8 1.0 



K/lc^n) 1 



FIG. 3: Excitation spectrum in the broken-axisymmetry 
phase. The solid curves represent the exact solutions and 
the dashed ones are approximate solutions given in Eq. (32). 
The linear and quadratic Zeeman coefficients are chosen to be 
p/|ci|n = 9/10 and q/\c\\n = 11/10. The modes (3 and 7 are 
the gapless modes associated with the simultaneous breakings 
of U(l) and SO(2) symmetries. The exact and approximate 
solutions for the /3 mode cannot be distinguished in this fig- 
ure. 



so that Eq. (18) reproduces Eq. (11). Therefore, we find 
that the normalization factors /i p , /i m , is p , and v m in 
Eq. (30), which determines & kj(T through Eq. (12), must 
be proportional to either e k ^ 4 or e\/ 4 . From numerical 

analysis, we find that they all have an e k 1//4 dependence 
as is the case in gapless excitation modes in the other 
phases (see Appendix B). It follows then from Eqs. (30) 
and (31) that (U + V)„, m ~ 0(e k 1/4 ) and (U - V) CT , m ~ 
0(e k /4 ), and 

(U + V) CT , m » (U - V) ff , m {tT = 0oi 7) (33) 

for e k — > 0. Therefore, we can neglect the second term in 
the square bracket in Eq. (17), obtaining 

k.a * ]T(U + V) ff , m (a Km + a f _ k , m ) (34) 

m 

for a = /3 and 7. The corresponding two Bogoliubov 
operators are then written as 



b kn ~ v v J2 m Cm (a k ,m + al Kr 



- Mm E m m Cm (ftk,m + al k , m ) , 



(35) 



r 



Equation (35) indicates that the quasiparti- 
cle operators for the two gapless modes con- 
sist of the number fluctuation operator SN — 



J2m Cm ( 



ation operator SF Z 



-k,. 



and the spin fluctu- 



J2 m m Cm (a k ,m + a 



-k,m 



We recall that the operator N is the generator of the 
global phase rotation and the operator F z is that of the 
spin rotation around the z axis. The creations of the 
gapless quasiparticles therefore lead to a change in the 
global phase of the order parameter and a rotation of the 
magnetization around the z axis. The modes f3 and 7 
are thus the Goldstonc modes that restore the U(l) and 
SO (2) symmetries. Since 8N and 5F Z can be regarded as 
the phonon and magnon operators, respectively, phonons 
and magnons are coupled in the quasiparticles described 
by Eq. (35). This is in contrast with the cases of the 
ferromagnetic and polar phases, in which phonons and 
magnons are decoupled (see Appendix B). The numeri- 
cally obtained coefficients [i p , [i m , and v m are shown 

in Fig. 4 as functions of q/\c\\n. We see that 6 k ,/3 is 
mostly the density fluctuation operator, while & kj7 is the 
linear combination of the number and spin fluctuation 
operators with roughly equal weights for small q/\c\\n. 
In other words, the /3-mode is a phonon-dominant mode 



and the 7-mode is a phonon-magnon coupled mode. The 
/3-mode crosses over to a phonon mode across the two 
neighboring phase boundaries, while the 7-mode crosses 
over to a magnon mode. 



D. Coherent excitations 



We investigate the dynamics of the states in which the 
quasiparticles are coherently excited. The excited state 
is assumed to be a coherent state 



(36) 



where |0)b is the vacuum of the Bogoliubov quasiparti- 
cles. The change in the expectation value of an observ- 
able Q due to the excitation of quasiparticles is given 
by 

(SQ^(t)) - WkMQn(t)Wk,,) -b<0|Qh(*)|0)b, (37) 

where the subscript H denotes the Heisenberg representa- 
tion. Since A k = T U(fc)B k — T V(fc)B_ k from the inverse 
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q/|ci|n q/|G|n 

Upfp-mode phonon) 

Vmiy-mode magnon) 

Vp (-/-mode phonon} 




(a) 



1.0 1.2 1.4 1.6 1.8 2.0 2.2 

q/|ci|n 



FIG. 4: (Color online) Coefficients in the phonon-magnon 
modes (fi p , v p , and ^ m in Eq. (35)) as functions of 

the normalized quadratic Zeeman energy q/\c\\n. The linear 
Zeeman energy is chosen to be p/\ci\n = 9/10. The verti- 
cal axis is scaled by (ek/|ci In) 1 / 4 /V No. The values of the 
scattering lengths obtained by Kempen et al. [20] are used: 
ao = 101.8 a.u. and = 100.4 a.u. The insets show the 
enlarged view of the curves for /x p and /x m - 



relation of Eq. (12), we obtain 



(** m (t))k,a 

_ \0k,* 



[(U - V)^™ cos(k • r - uj a t + 4> k ^) 

+ i(U + V) a ,m sin(k • r - + 0k jO -)], 

(38) 



where 4>k.a = arg(/3fc j(T ) and u> a = E^ a /h. Since the ra- 
tio of the real part to the imaginary part is estimated 
from Eq. (33) to be (ek) 1 ^ 2 <C 1, the real part is neg- 
ligible for the two gapless modes, a — [3 and 7, in the 
long- wavelength limit. Therefore, {6^ m (t))k >(T is almost 
entirely imaginary, which indicates that the change oc- 
curs mostly in the phase of the real order parameter Q m . 
Thus, the excitations of (3 and 7 modes lead to a global 
phase rotation and a spin rotation around the z axis. 

To study how the quasiparticle excitation rotates the 
spin, we calculate (<5F(i))k, CT - Keeping terms up to those 
of the first order in (3k m we obtain 



V 1 * 



'*>/V±H.c. 



(39) 



where (£ = x. y, z) are the spin-1 matrices defined 
in Eq. (A13), and the upper signs refer to £ = x and 
z and the lower signs to £ = y. Hence, Eq. (33) leads 
to (5F y )^ a > (SF x ) Kcr and (5F v ) k ,* > (5F Z ) K *- This 
is because (F x ) ^ 0, (F z ) ^ 0, and (F y ) = from the 
assumption of real and positive Q m , and the infinitesimal 
spin rotation around the z axis changes only (F y ). 
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FIG. 5: Oscillations of the order-parameter phases (xm = 
tan -1 [Im(-(/) m )/Re('i/) m )]) and the azimuthal angle of magne- 
tization around the z axis (<j> = tan -1 [(F V )/{F X )]) caused by 
the excitation of the Bogoliubov quasiparticles in the long- 
wavelength limit. The vertical axis is marked in arbitrary 
units. The Zeeman energies are assumed to bep/|ci|n = 9/10 
and g/|ci|n = 11/10. (a) Phonon-dominant mode (/3-mode). 
Since rotations of the order-parameter phases are much larger 
than that of the magnetization, we cannot see the latter in 
this figure, (b) Phonon-magnon coupled mode (7-mode). 



Thus, we have shown that the excitations of the Gold- 
stone modes (3 and 7 lead to U(l) and SO(2) transfor- 
mations. Oscillations of the order-parameter phases and 
those of the azimuthal angle of the spin vector are shown 
in Fig. 5. Figure 5(a) shows that the /3-mode excitation 

changes the phases of (V>i), (ipo), an d (^-1) in the same 
manner. This is because as shown in Fig. 4 the dominant 
contribution to the (3 mode is made by phonons which are 
insensitive to individual spin components. On the other 
hand, the 7-mode excitation describes not only fluctua- 
tions of the overall phase but also those of the spin vector 
around the z axis. Since the rotation around the z axis 
is e l ' pFz \m) = e lipm \m), xi and X-i are out of phase with 
respect to xo- 

The gapful mode (a-mode) can be interpreted as play- 
ing the role of changing the magnitude of magnetization. 
As shown in Fig. 6, the fluctuation of F x is dominant 
when the a-mode is excited. The z component F z can- 
not vary due to the spin conservation, and hence the spin 
fluctuation is restricted in the x-y plane as illustrated in 
the inset of Fig. 6. 



IV. DISCUSSION 

We have shown that the ground state in the shaded 
region of Fig. 1 is the broken-axisymmetry phase featur- 
ing transverse magnetization. Here, we discuss possible 
experimental consequences of the axisymmetry breaking 
and the transverse magnetization. 
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FIG. 6: Oscillation of transverse magnetization (F x ) due to 
coherent excitation of the gapful mode (a-mode). The inset 
schematically illustrates the change in the spin vector caused 
by the excitation of the a-mode. 



Let us consider a situation in which atoms are pre- 
pared in the m = state. Whenp = and q < 2|ci|ra, the 
ground state is in the broken-axisymmetry phase and the 
transverse magnetization should develop in time. How- 
ever, the total spin angular momentum parallel to the 
magnetic field must be conserved, and for small q the 
magnitude of the transverse component of the total spin 
is nearly entirely conserved. Consequently, local magne- 
tization varies in space only insofar as the total spin is 
conserved [23, 24]. This constraint leads to formation of 
various spin textures depending on the trap geometry. 
For example, in an elongated cigar-shaped trap, a stag- 
gered domain structure or a helical structure is formed 
spontaneously [23]. In a pancake-shaped trap, on the 
other hand, a concentric domain structure is formed [23]. 

The formation of the domain structure costs the kinetic 
energy and the ferromagnetic energy at the domain walls. 
If the direction of the spin vector changes gradually over 
space, the formation of the spin texture costs little en- 
ergy. One of such textures is a topological spin texture, 
in which the orientation of transverse magnetization has 
a 27r winding about a central defect. When the size of 
the system is small, the topological spin texture becomes 
energetically favorable and develops spontaneously [24]. 
Recently, formation of the topological spin texture has 
been observed by the Berkeley group [25], in which the 
state of the system is changed rapidly by a change in the 
magnetic field from the |0) region to the shaded region in 
Fig. 1. The spontaneous transverse magnetization in the 
Berkeley experiment is a manifestation of the symmetry 
breaking discussed in the present paper. If the amount of 
the change in the magnetic field is small or if the speed of 
the change is slow, the system is not markedly disturbed, 
and low-energy gapless excitations (the (3 and 7 modes) 
should be observed. 



V. CONCLUSIONS 

We have studied a spin-1 ferromagnetic BEC by tak- 
ing the quadratic Zeeman effect into account. The mean 



field theory predicts that BECs with ferromagnetic in- 
teractions show the broken-axisymmetry phase, in which 
magnetization tilts against the direction of an external 
magnetic field. Here, the SO (2) symmetry is broken, in 
addition to the U(l) global phase symmetry. Applying 
the Bogoliubov theory for a BEC with spin degrees of 
freedom, we have found one gapful mode and two gap- 
less Goldstone modes for this phase. We have analytically 
shown that two gapless modes are the coupled phonon- 
magnon modes that restore the U(l) and SO(2) sym- 
metries simultaneously. Numerical analysis has shown 
that one Goldstone mode is the phonon-dominant mode 
and the other is the phonon-magnon coupled mode with 
roughly equal weights. The gapful mode changes the 
length of the spin by fluctuating the spin in the direc- 
tion perpendicular to the magnetic field (see the inset of 
Fig. 6). 

When more than one continuous symmetry is sponta- 
neously broken, multiple gapless modes emerge and they 
couple with each other to form Goldstone modes as shown 
in this paper. Such multiple Goldstone modes may there- 
fore be found in spin-2 BECs [12] and higher spin BECs, 
which merit further investigation. 
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APPENDIX A: DERIVATION OF THE PHASE 
DIAGRAM FOR ci < 

In this appendix, we derive the analytic expression 
for the spin components in Eq. (6) in the broken- 
axisymmetry phase and reproduce the phase diagram for 
ci < as shown in Fig. 1. 

The average energy per atom is given from Eq. (5) by, 

( Wbec ) 
e = — 

N 

= ^|2CiC-i-C 2 | 2 + E (-pm + qm 2 )\U\ 2 . 

(Al) 

The first term on the right-hand side of Eq. (Al) can be 
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rewritten as 



— — I^CiC— i - Co | 
kiln 



2|CiC-i| - ICol'e 



2„i(2xo-Xi-X-i) 



(A2) 



where Cm = |Cml eiXm - Hence, the phase depends on en- 
ergy only through 2xq — Xi ~ X— 1> and the energy is mini- 
mized for 2\o ~ Xi — X-i = 0. Without loss of generality 
we assume that Cm's are real and non-negative. 
We minimize 



K = e- + $ + 



(A3) 



where a Lagrange multiplier /j, is introduced to ensure the 
normalization condition Ci + Co + C-i = 1- Stationary 
conditions can be obtained through differentiation of K 
with respect to Cm' s as 



dK 

9c±T 

dK 
5Co" 



= 2| Cl |n(2CiC-i-Co 2 )CTi + 2(gTp-M)C±i =0, 



= -2| Cl |n(2CiC-i-Co 2 )Co-2Ko = 0. 



ft follows from Eq. (A5) that either Co = or 
Ai =-|c 1 |n(2CiC-i-Co 2 ) 



(A4) 
(A5) 

(A6) 



has to be satisfied. When Co = 0, we find that K is 
minimized with 



Ci = l, C-i=0 forp>0 {h = q-p), 
Ci=0, C-i = 1 forp<0 (h = q + p). 

When Co 7^ 0, Eq. (A4) becomes 

(?TP-/*)C±i =K=fi- 
We can easily see that the polar state 

Co = 1, Ci = C-i = (h = i) 



(A7) 
(A8) 



(A9) 



(A10) 



satisfies Eq. (A9). Using Eqs. (A6) and (A9), we obtain 
the solution for the broken-axisymmctry phase in Eq. (6). 
This solution is valid for 



Q 



p 2 > and 2\ Cl \nq - q 2 + p 2 > 0. 



(All) 



The energy of the broken-axisymmetry state is calculated 
to be 



ebr 



A{\ci\nq)' 



;(q 2 -p 2 )(p 2 -4\ Cl \nq-q 2 ). (A12) 



We can show that e br < 1 and et, r < q ± P are always 
satisfied in the region of Eq. (All), and thus we obtain 
the phase diagram in Fig. 1. 



by 



The three components of the spin-1 matrices are given 




fz = I I , 

0-1 



(A13) 



and the corresponding spin components are given by 

(f x ) = No^uc = ^ [Co (Ci + c-i) + Co (cr + CO] . 

(AW) 

(F y ) = N tff y C = ^ [Co (Ci - C-i) - Co (Ci* - C-i)] , 



y/2 

(F z ) ee AoCt/.C-^odCil'-IC-il 2 )- 



(A15) 
(A16) 



Substituting Eq. (6) into Eqs. (A14)-(A16), we obtain 
Eqs. (7)-(9). 



APPENDIX B: EXCITATION SPECTRA IN THE 
FERROMAGNETIC AND POLAR PHASES 

1. Ferromagnetic phase 

We first consider the ferromagnetic phase (| + 1) in 
Fig. 1) for p > 0. In this phase the order parameters are 
given by Ci = 1 and Co = C-i = 0. The matrices M and 
N in Eq. (19) are shown to be 

' e k + cqji + c\n 
M = I e k • P q 

e k + 2p - 2cin , 



(Bl) 
(B2) 



cqh + c\n 0^ 
X = I 
0, 

and the matrix G in Eq. (22) is given by 



Ak(ek + 2g 2 n) 

G= (e k +p-q) 2 

\ (e k - 2p - 2c x n) 2 

(B3) 

Diagonalizing G we obtain three excitation energies as 

E p = v /e k (e k + 2g 2 n), 
E = e k +p-q, (B4) 
E-x = e k - 2p-2cm, 
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and the associated quasiparticle operators as 



6k >p = -rm a ^ + jm a - 




k,l' 



(B5) 



1/2 

portional to e k and 6k, p in Eq. (B5) is proportional to 
e k 1 ^ 4- The singular ek-dependence in each of .E p and 6k, p 
is therefore canceled in the product Epb^ p 6 k>p , giving the 
original ek-dependcncc in the effective Hamiltonian. 



The corresponding results for p < can be obtained by 
changing the linear Zeeman energy p to —p. For q = 
these results reduce to those obtained in the absence of 
the quadratic Zeeman effect [17]. The excitation spec- 
trum E p is similar to the one obtained for scalar BECs; 
it is a phonon mode with the speed of sound given by 
c = y / g 2 n/M. The other two modes are understood 
as excitations from the m = 1 state to the m — and 
m = — 1 states, respectively. In the long- wavelength limit 
(k — > 0), Eq (E_i) coincides with the single-particle en- 
ergy difference between the m = 1 state and the m = 
(m = — 1) state. Hence, the phonon mode (E p ) and the 
magnon modes (E ,E_i) are decoupled in the ferromag- 
netic phase. 

It can be easily checked that these excitation modes 
and operators do reproduce the effective Hamiltonian 
(11) through Eq. (18). For small ek, E p in Eq. (B4) is pro- 



2. Polar phase 

In the polar phase (|0) in Fig. 1), the order parameters 
are Co = 1 and Ci = C-i =0, and the matrices M and N 
have the forms 

f £k -p + q + c\n 
\J i . k • ':,n 

ew+P + q + cxn, 



cm^ 

S = i c n 
iCin 

which give 



(B6) 
(B7) 



'( £ k — P + q + c\n) 2 — c\n 2 2c\np 

G = I e k (e k + 2c n) 

— 1c\np (ek + p + q + c\n) 2 — c\n 2 



(B8) 



The eigenvalues are obtained as 



E p = v/ £ k(ek + 2c n), 

= -p + 6c, 



(B9) 



where £ k = \J (ek + <?)( e k + 9 + 2cin). The quasiparticle operators are found to be 



fok, P = 



y^k + c n + Ep 



ak,o + 



\Ak + can - Ep 



t 



-k,0' 



V^k + 9 + cm + $ k . Vek + g + cin-^ k . t 
0k,± = 7=== a k,±l === a 



k,Tl- 



(BIO) 



The first mode £p in Eq. (B9) is a phonon mode. The Thus, also in this phase, phonons and magnons are de- 
other two are magnon modes which scatter a quasiparti- coupled, 
cle from the state with m — to the one with m = ±1. 
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